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We propose a LHC search for dilepton resonances in association with large missing 
energy as a generic probe of TeV dark sector models. Such resonances can occur 
if the dark sector includes a 17(1) gauge boson, or Z', which kinetically mixes with 
the Standard Model 17(1). For small mixing, direct 2—7-1 production of the Z' is 
not visible in standard resonance searches due to the large Drell-Yan background. 
However, there may be significant production of the Z' boson in processes involving 
other dark sector particles, resulting in final states with a Z' resonance and missing 
transverse momentum. Examples of such processes include cascade decays within the 
dark sector and radiation of the Z' off of final state dark sector particles. Even when 
the rate to produce a Z' boson in a dark sector process is suppressed, this channel 
can provide better sensitivity than traditional collider probes of dark sectors such as 
monojet searches. We find that data from the 8 TeV LHC run can be interpreted to 
give bounds on such processes; more optimized searches could extend the sensitivity 
and continue to probe these models in the Run H data. 
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I. INTRODUCTION 

There is compelling evidence that most of the matter in the Universe is composed of 
nonbaryonic particles, the dark matter (DM), the nature of which is otherwise unknown. 
One of the only quantitative data points known about dark matter is its current cosmological 
energy density, 0 _dm ~ 0.25 pQ. If the dark matter is a relic of freeze-out from the thermal 
bath of the early Universe, then this relic density can be achieved |2] with a dark matter 
mass of 0(100 GeV) with couplings of order ~ 0.1. This “WIMP miracle” is a prime 
motivation for considering dark matter production at the LHC. Dark matter candidates 
of weak scale mass can also naturally emerge within theories developed to address the 
electroweak hierarchy problem, such as supersymmetric models. 

Current searches for dark matter at the LHC focus on pair production of invisible DM 
particles plus radiation from the initial state in the form of jets, photons, or electroweak 
bosons [3H22]. The resulting “monojet”, “monophoton” etc. signatures have considerable 
SM background, but still allow for constraints to be placed on dark matter interactions with 
quarks and gluons. Within an effective held theory (EFT) framework, these processes can 
be correlated with signals from elastic DM scattering in detectors (direct detection) and 
astrophysical DM annihilation (indirect detection) [23H2H]- 

This program is appropriate for the minimal assumption of a single DM particle and no 
other new physics. However, when one goes beyond this minimal framework, other types 
of collider searches may provide much more powerful probes [2SH3H]- A familiar example 
is that of supersymmetric models, in which searching for squarks and gluinos decaying to 
a neutralino DM candidate is usually a far more effective probe of the new physics than 
searches for direct neutralino production. More generally, any new particles associated 
with dark matter can provide additional collider signatures which may greatly enhance the 
prospects for discovery. Note that the same enhancement does not extend to direct and 
indirect detection of dark matter, which are only sensitive to the actual cosmological relics. 
In this respect colliders provide a unique window into the physics associated with dark 
matter. 

In this work we specialize to the case where all new particles are gauge singlets under 
the Standard Model, forming a “dark sector.” Although the model space for such a dark 
sector is vast, some well-motivated assumptions greatly narrow down the possible collider 
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phenomenology. Consistent with renormalizable held theory, we can consider the new par¬ 
ticles to be either fermions, scalars or gauge bosons. Since the coupling of these states to 
the SM is generally weak, once dark sector particles are produced at colliders they will tend 
to cascade decay within the dark sector until that is no longer kinematically possible. In 
particular, particles with arbitrarily weak couplings to the SM can still be produced through 
decay of or radiation off of other dark sector particles. If the theory preserves baryon and 
lepton number, or an analogous dark fermion number, then the lightest dark fermion will be 
absolutely stable and appear as missing energy at colliders. Dark bosons however may not 
be protected by any quantum numbers and could decay into the Standard Model. Although 
this decay width may be small due to weak couplings, the branching ratio of the dark boson 
to the SM can still be large if decays to other states are suppressed- particularly if they 
are not kinematically accessible. This allows for visible collider signatures from extremely 
weakly coupled particles, as in “hidden valley” models [lOl ITT] . 

The collider signatures of a dark gauge boson, or Z\ can be particularly striking. The 
most relevant possible interaction of a dark vector boson with the Standard Model is through 
the kinetic mixing “portal”, i.e. the operator where is the SM hypercharge held 

strength and is the dark gauge boson held strength |12]. As a result of this mixing the 
Z' will couple to SM currents and thereby decay into pairs of fermions, including an 0(1) 
branching ratio to leptons. This gives a distinctive and easily measured dilepton resonance 
in collider events. Direct 2 —)■ 1 production of Z's has been searched for extensively at 
colliders [^3] - H6] , but the SM Drell-Yan background limits the sensitivity of such searches to 
weakly coupled Z’s. However, dark sector cascades will tend to produce Z's in association 
with missing transverse momentum (MET) from invisible dark fermions. In this work we 
show that leveraging this large MET would allow for searches with very low background 
and high acceptance for a broad class of dark sector models. We give examples of simplihed 
models, representative of more general possible dark sectors, which could be discovered by 
such an analysis where current approaches such as mono jet searches and inclusive resonance 
searches fail. 

In the next section we describe a proposed search for dilepton resonances plus MET, 


building off of existing analyses by ATLAS and CMS. In section we discuss a variety 
of signal models (see e.g. hgure which realize the resonance plus MET hnal state and 
demonstrate the reach of the proposed search for each of them. 
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FIG. 1: Feynman diagrams for some of the signal processes we discuss in section 


III 


In the 


process of £g. la, which we refer to as “darkstrahlung” (section III A), a Z' boson is 


radiated off of hnal state dark sector particles. In the Z' is produced in the cascade 
decay of one dark sector particle to another (section |IIIB ). For both of these processes, we 
will consider production of the initial dark sector particles through a heavy mediator which 


can be integrated out. In £g. Ic, an on-shell scalar 0 is produced through mixing with the 


SM Higgs, and decays to two Z' bosons (section IIIC). 


II. DILEPTONS PLUS MET AT THE LHC 


A search for dilepton resonances in association with missing transverse momentum (MET) 
can be understood as a straightforward extension of existing, non-resonant LHC studies of 
dileptons plus MET. Within the SM, analyses in this hnal states are used to measure the 
cross section of W~^W~ ■ UZHIH]. Searches in this channel can also probe new physics such 
as decays of charginos and sleptons in supersymmetric models [50]. In this work, we use 
cuts building off of the 8 TeV ATLAS search for charginos [Sn|, in particular their signal 
region WWa. The signal region requires two same havor, opposite charge leptons. The cuts 
for leptons, missing energy, and jets for the signal region are given in Table (all objects 
are defined after subtraction of tracks and calorimeter energy from pileup interactions, and 
standard isolation cuts are imposed for leptons). The missing energy variable, J^Trei^ 
dehned as 
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Cuts 

PT,e 

> 10 GeV 

\Ve\ 

< 2.47 

PT,fi 

> 10 GeV 

I^mI 

< 2.4 

PT,leading lepton 

> 35 GeV 

PT,second lepton 

> 20 GeV 

'mu' 

> 20 GeV, < 120 GeV 

\mu' - mz\ 

> 10 GeV 

Pt,u' 

> 80 GeV 

^Trel 

> 80 GeV 

jet veto 

events with pj > 20 GeV and \r]j\ < 2.4 
events with pJ > 30 GeV and 2.4 < \r]j\ < 4.5 


TABLE I: The cuts employed by ATLAS for signal region WWa in their search for SUSY 
electroweakinos in 20 fb~^ of 8 TeV data [50]. The analysis requires exactly two opposite 
sign leptons. The signal processes are shown in Fig. while the main backgrounds are the 
leptonic decays of and ti. Inside the Z mass window, ZZ process also contributes 

to the background. In the table above, mz is the mass of SM Z-boson. I^Trei is defined in 
the text. Our proposed searches for Z' + MET signals involve changes to the cuts shown 

in bold. 


^T,rel - 


( 1 ) 


l^T if ^(j^i > 71 /2 

]^TsinA(j)i if A0£ < 7r/2 
where Acpg is the azimuthal angle between the I^t with the nearest lepton or jets. 

Although the ATLAS signal selection imposes only the cuts on mu listed in Table 13 full 
distributions of the expected and observed values after all other cuts are shown in ref. 


(see Figure 3a of [50] or Fig. 12 of this work). From this information one can derive bounds 


on resonant dilepton plus MET production for a range of masses. However, the mu data 
as presented is binned in intervals of 10 GeV, which is not optimal for searching for narrow 
(weakly coupled) resonances. The lepton energy resolution at ATLAS ranges between 1% 
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to 4% depending on the lepton flavor and rapidity [5T1 [52]. For CMS, the energy resolution 
ranges between 1% to 6% [SnilMI- Given these resolutions, we expect that the dilepton mass 
resolution to be better than < 3% (52]. As a comparison, an ATLAS simulation of Higgs 
decaying to a dimuon pair [55] estimates a FWHM of 5.6 GeV (la resolution of 2.4%) for 
\r]fj_\ < 1. For the purpose of this paper, we choose the mass window of mz' ±2.5%m^/ when 
estimating the current optimized and future bounds. 

Depending on the signal model, altering the cuts on pT,ei' and l^Trei could also improve 
the sensitivity. As we will show in the next section, certain models will tend to produce 
Z'’s with large boosts. For such models, higher values for the pr/e and l^Trei cuts are more 
desirable, as one can eliminate backgrounds almost completely without having a significant 
reduction in signal efficiency. In other scenarios, however, the Z' can produced with low 
velocity, so that a looser cut on pr/t and l^Trei is then necessary to have appreciable signal 
acceptance. Therefore in a broad search for new resonances in association with missing 
energy, there is motivation for multiple signal regions with varying cuts. For simplicity, in 
this work we will always take the same numerical value for the Pt,u> cut and the l^Trei cut, 
and refer to this value as simply the “MET cut.” 


III. SIGNAL MODELS 


Before proceeding to specific models of Z' + MET production at the LHG, we review the 
basic aspects of a kinetically mixed Z'. We consider a dark gauge boson of a group U{1)d, 
which couples to the SM U{1)y via the kinetic mixing operator. The Lagrangian of the 
model is given by 

+ + ( 2 ) 

where is the dark gauge boson, is the field strength tensor of X^, and is the 
field strength tensor G(l)y. These gauge bosons are written in their gauge eigenstates. 
ruz' is the mass of dark gauge boson, which can arise from various mechanisms, such as 
spontaneous breaking by a dark Higgs or the Stueckelberg mechanism. The mixing term 
can be generated by a loop of heavy particles charged under both U (l)y and U{1)d- In this 
scenario, an estimate of the kinetic mixing e is given by [56] 


e ~ 



(3) 
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where go is the U{1)d coupling constant, m is the mass of the heavy particle that is charged 
under both gauge groups and A is a cutoff scale. A wide range of mixing values can be 
realized in various scenarios: e 10-12 tQ ^ _ iq-3 

Existing experimental data place various constraints on Z' models. A summary of the 
bounds for a wide range of masses and kinetic mixings can be found in [59l |66] . However, 
most of the bounds derived in these references apply only in the case where the Z' decays 
dominantly to the visible sector. In a model-independent bound is derived by con¬ 
sidering the effect of the Z' mixing on precision Z boson measurements at LEP [HHl EH]- 
For the range of Z' masses that we are concerned with, ~ 10 — 1000 GeV, the model- 
independent bounds on the mixing parameter are typically e < 0.01 — 0.1 (green region in 
figure . 

Given this bound on the mixing, there is still considerable parameter space in which the 
Z' can decay promptly to the SM. When the Z' decays only to the SM sector, the branching 
fraction to each decay channel is independent of e and other parameters in the dark sector, 
and depends only on the Z' mass (for small e). In Figure]^ we plot the branching fraction 
of the Z' to leptons (e+e” and in this scenaric0 

Another possibility is that the Z' can decay to dark sector states. The branching fraction 
in this case does depend on the kinetic mixing e and also the details of the dark sector, i.e. 
the dark gauge coupling constant and the spectrum of states charged under the dark f/(l). 


In Fig. 2b, we show contours of the branching fraction of Z' —)■ in the case of only 

one light dark sector fermion y that the Z' can decay to. As a benchmark point we take 
niy. = 10 GeV and go = 0.1 (and unit charge for the y). 


A. Darkstrahlung 


One generic way to produce a Z' in a collider event is to radiate it off of hnal state 


dark sector particles (hgure la), a process which we will refer to as “darkstrahlung.” To 


illustrate this, we consider a simplified model of a dark sector fermion y that is produced 
^ The decay width of Z' to leptons, ignoring the lepton mass, is given by T = 
^ (^(^gcwSa - 35' (^swSa “ + {q So. +9' (^swSa “ ^ > where g and g' are the 

SU(2)w and U(l)v coupling constants respectively, sw and cw are the sine and cosine of the weak 
mixing angle respectively. Sq, and Cq, are the sine and cosine of the mixing angle between the Z', Z and 
7. In the limit of < 1, is given by . 









Branching fraction of Z rr for = ^0 GeV and g© = 0.1 



50 100 150 200 250 300 350 400 

[GeV] 


(a) > raz' 



FIG. 2: Left: Branching fraction of the Z’ to two leptons (e’''e“ and as a fnnction of 

mass, assnming the Z' decays only to the SM sector. Right: Contonrs of the same 
branching ratio, in the case where there is a 10 GeV dark fermion which the Z' can decay 
with conpling strength gjj = 0.1. The green region on the hgnre is exclnded by precision Z 

observables at LEP |53- 

at hadron colliders throngh an effective operator of the form ^ {Q'lfj.'lbO) (X7^75X)0 When 
initial state radiation is taken into acconnt, snch models can be probed by searches for 
monojets, monophotons, etc. (SHH]- However, if the x particles are charged nnder a U{1) 
gange symmetry, then they too can radiate, prodncing a Z' boson in the hnal state, with 
probability determined by the dark U{1) conpling. As we have discnssed, snch a Z' will 
always decay to Standard Model states if the dark sector is not kinematically accessible 
(here, if mz' < leading to an observable hnal state resonance in association with large 

missing energy from the stable x particles. 

The darkstrahlnng process can easily be the dominant prodnction mode of the Z' at 
colliders. Direct 2 —)■ 1 prodnction of the Z' is snppressed by the kinetic mixing parameter 
e, and for e ;< 10“^ [TOl [7T] is too small to be discover amidst the Drell-Yan backgronnd. In 
^ While the exact Lorentz structure of this operator has little effect on collider phenomenology, if x is a cos¬ 
mological relic then its non-relativistic scattering through this operator generates mostly spin-dependent 
interactions in direct detection experiments, which are relatively unconstrained. 
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FIG. 3: The production cross section for pp —>■ XX^' the 8 TeV LHC as a function of 
for = 200 GeV, A = 1 TeV and = 0.1. (These parameters are just out of reach 
of current monojet searches, see Figure and associated discussion.) 

contrast, the cross-section for the darkstrahlung process depends only on the magnitude of 
the qqxx operator and the dark U{1) coupling, and the branching ratio for Z' —)■ £~^£~ is 
independent of e if dark sector decay modes are inaccessible. Furthermore, the missing energy 
signature in this channel allows one to almost completely reject Drell-Yan background, 
greatly increasing the signal sensitivity. 

Because the dark matter production proceeds through an operator suppressed by a high 
scale A, most events occur at high center-of-mass energy ^/s, resulting in the y particles 
having large boosts. The high ^/s enhances the energy of hnal state radiation as well, leading 
to high transverse momenta for the Z'. The rate of FSR includes an enhancement by s/m\, 
due to the longitudinal mode of the Z' (similar behavior is seen e.g. in bremsstrahlung of 
electroweak bosons in the annihilation of heavy dark matter particles [ZH [73]). Figure 
shows the darkstrahlung cross-section as a function of Z' mass for an example set of param¬ 
eters. Figure]^ shows the MET distribution for darkstrahlung events, compared to Standard 
Model WW events with the same dilepton invariant mass. 

The kinematics of darkstrahlung events depend only on the masses of the y and Z' 
particles; we can therefore interpret collider searches as placing bounds on the darkstrahlung 
cross-section as a function of these two parameters. As discussed in section [TT| results from 
the ATLAS search for electroweakinos can be used to place some constraints on this model. 
Using the data from Figure 3a of [501, which is given in dilepton invariant mass bins of 10 
GeV, we derive constraints on the darkstrahlung cross-section times branching ratio to £'^£~ 
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g^missing 

FIG. 4: Comparison of the MET distributions for the darkstrahlung process with 
= 200 GeV, = 50 GeV and the WW —)■ t^l~vv background with the dilepton 
invariant mass hxed to mu' = 50 ± 2.5 GeV. The cuts for both darkstrahlung and WW 
processes are the same with the ATLAS cuts shown in Table [T| except we require 
\mu' — 50 GeV| < 2.5 GeV and the J^Trei not used. 


shown in hgure 5a (Details of our MG simulation can be found in appendix]^) A more 
optimized analysis would use smaller invariant mass bins as well as a slightly higher MET 


cut to further reduce the SM background. In hgure we show the projected reach for an 
improved search making use of these tighter cuts, giving an 0(1) improvement. 

In terms of the absolute cross-sections probed, these analyses are far more sensitive to 
this model than monojet searches, which constrain cross-sections of 0(10 pb). However, the 
rate for darkstrahlung depends on the dark U{1) coupling and at weak coupling can be much 
lower than the monojet rate, though even then the darkstrahlung process can be a better 
probe due to the much greater sensitivity. For large dark 0(1) couplings, there can be an 
appreciable rate to radiate two or more Z'^s in a dark sector event, in which case four lepton 
-|- MET signatures can be produced. Such hnal states have close to zero SM background 
even without reconstructing the Z' resonances [71] and can also be used to place constraints. 
Figure]^ shows how the discovery reach of monojet searches, our dilepton -|- MET proposed 
search, and multilepton searches compare as a function of the dark U{1) coupling For these 
parameters, a dilepton -|- MET search outperforms mono jet searches even for as small as 
~ few X 10“*^, and is superior to a multilepton search even at strong coupling od ~ 1- There 
^ The ZAMET bound shown in Figure]^ can be approximated as < {6.8/an) fb, while the multi- 

pp^xx ^ (23/o|)) fb. 


lepton bound is given by cTj 
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ATLAS bounds on (crx BR)pp^ /+ 
(Darkstrahlung Model) 



5 


100 200 300 400 500 

(GeV) 

(a) 


Projected optimized reach in (crx BR)pp ^ /+ /- ^^(fb) 


(Darkstrahlung Model) 



100 200 300 400 500 


(GeV) 

(b) 


FIG. 5: Left: 95% C.L. constraints on the cross-section times branching ratio of 
PP XX^' 1 for the darkstrahlnng model in the plane of x mass and Z' mass, 

from data presented in an ATLAS search for leptons and missing energy (Figure 3a of 
ref. ISO]). In the gray shaded region, the Z' is able to decay into x particles instead of 
being forced to decay to SM states. Right: Projected 95% exclusion reach of a more 
optimized search, with the MET cut increased to 100 GeV and the dilepton invariant mass 

window narrowed to ±2.5% of the Z' mass. 


is considerable parameter space for which the production operator ^ (^ 7 ^ 750 ') (X7^75X) can 
be probed to A as high as a few TeV. For these values of A, EFT validity at LHG energies is 
achieved even if the UV theory has order one dimensionless couplings, unlike the parameter 
space probed by monojet searches [7SHZZ!. 

The model we have discussed is presented merely to illustrate possible features of dark 
sector collider events. A complete model of dark matter would have to address other issues 
such as the relic abundance. The introduction of a dark (7(1) offers a way to produce an 
appropriate thermal relic abundance; if > mz' then the x particles can annihilate to 
Z'’s which then decay to the Standard Model. This can give the correct DM abundance 
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ao 

FIG. 6: Comparison between bounds from monojet [3], Z^'+MET and multi lepton m as 
a function of the dark gauge coupling = 5'D/(4vr) for = 200 GeV, mz' = 120 GeV. 

Multilepton bound is always worse than Z'+MET bound for «£> < 1. 

given weak scale masses in the dark sector and couplings of 0(0.1) (e.g. = 200 GeV, 

mz' = 120 GeV, «£> = 6 x 10“^), realizing the “WIMP miracle.” If the y is identified as 
dark matter then it can give direct detection signals mediated by the same EFT operator 
we used to describe collider production. For Dirac dark matter the vector-channel operator 
(qIiiQ) ixi^x) is also present and gives large spin-independent direct detection rates, but 
for Majorana fermions this operator is absent. (In the Majorana case the dark matter mass 
term can then be generated from a Higgs which breaks the dark f/(l).) Of course, all of the 
results regarding collider physics hold true even if x is not the cosmological dark matter, 
e.g. if it decays to other dark sector particles, if it decays to the SM with a long lifetime, 
or if it has a small relic abundance. This reflects the fact that collider events can probe the 
dark sector even beyond dark matter. 


B. Cascade decays 

Another potential production mode for Z' bosons is in the decay of one dark sector state 
to another. A generic dark sector process may involve a number of such decays, with various 
branching ratios to hnal states, analogous to the rich set of hnal states from e.g. top or Higgs 
production in the Standard Model. As we argued in the introduction, dark sector cascades 
will typically produce large missing energy from some number of stable invisible particles. 
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but could also result in one or more bosonic resonances, such as a Z', which can decay 
promptly to the SM. 

Once again, a simple toy model can illustrate the generic phenomenology of this scenario. 
Consider extending the spectrum of the previous section with an additional dark sector 
fermion, X 2 , which can decay to a Z' and the lighter fermion, which we denote xi this 
section]^ Production of X 1 X 2 or X 2 X 2 can then lead to states with one or more Z'’s in 
association with missing energy, with rates depending on the branching ratio of X 2 Z'xi- 
As before, the Z' may be forced to decay to the SM if it is too light to decay to Xi- Final states 
with two or more Z' bosons are best probed by existing searches in multilepton channels, 
as discussed above, so we will focus on the hnal state with one Z'. Final states with one Z' 
can dominate if the branching ratio of X 2 —t Z'xi is small, or if X 1 X 2 production dominates 
over X 2 X 2 - Both the “off-diagonal” production X 1 X 2 and the decay X 2 —t Z'xi can naturally 
arise from spontaneous gauge symmetry breaking in models with Dirac fermions or complex 
scalars [7H] . 

For the purposes of our Monte Carlo study, we specialize to the case of X 1 X 2 production 
through the operator ^ (? 7 /i 75 <?) ( 717 ^ 75 X 2 )• Figure shows the expected reach in cross- 
section times branching ratio of the optimized 8 TeV search discussed in the previous section, 
with the Z' mass hxed to 40 GeV while the xi and X 2 masses are varied. The reach is 
approximately constant over the entire parameter space shown, even in the “squeezed limit” 
where ^ +mz', where the Z' is produced nearly at rest in the frame of the X 2 - The 

invariance of the collider signals with respect to the y masses is a result of the large boost 
of the y particles in this production model, which ensures that the Z' from the decay has 
high pt independent of the energy released in the decay itself. The cross section bound of 
~ 2 fb corresponds to a value of A of about ~ 7 — 9 TeV, depending on the mass spectrum. 

C. Dark Higgs 

Thus far we have not discussed the “dark Higgs” held <F that must be introduced to 
spontaneously break the dark 17(1) gauge symmetry in a weakly coupled model. Such a 
held can couple to the Standard Model at the renormalizable level through the “Higgs por- 


Although the Standard Model famously lacks such “flavor-changing neutral currents” at tree level, they 
can occur in more general theories, such as those with vectorlike fermions. 
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Reach in (crx BR)pp ^ /+ /- 
(Cascade Model, mz'= 40 GeV) 
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60 80 100 120 140 160 180 200 

m;^2(GeV) 

FIG. 7: Projected 95% exclusion reach in the plane of Xi and X 2 masses for the cascade 
decay model with = 40 GeV, from an optimized 8 TeV search as in hgure In the 
red region X 2 —t Z'xi is kinematically forbidden, while in the gray shaded region the Z' is 
able to decay into y particles instead of being forced to decay to SM states. 

tal” operator, |i7p|<l>p. When H and $ each acquire their vacuum expectation values, this 
operator induces mixing between the Standard Model and dark sector Higgses. This allows 
additional processes to create dark sector states at colliders, with very different kinematics 
compared to production through the higher-dimension operators discussed above. Never¬ 
theless, we hnd that searching for resonances accompanied by missing energy can once again 
probe much of the parameter space of generic models. 

We will focus on the production of a single on-shell scalar, denoted as 0, which mixes 
with the SM Higgs. This state can then be produced through all of the same channels as the 
SM Higgs; for simplicity we will focus on the dominant gluon-fusion-initiated 2 —1 process. 
Once produced, 0 may be able to decay to two Z' bosons, depending on the parameters of 
the model. If these Z' states are kinematically forced to decay to the Standard Model (as we 
assumed in the previous subsection), then there can be spectacular four-lepton hnal states 
reconstructing multiple resonances, as considered in [79II83] . However, if there are dark 
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8 TeV LHC, detector level 



FIG. 8: Comparison of the MET distributions for the (j) —)■ Z'Z' —)■ process with 

mz' = 50 GeV and two different values for the 0 mass, and for the WW —)■ 
background with the dilepton invariant mass hxed to = 50 ± 2.5 GeV. The cuts for 
both signal and WW processes are the same with the ATLAS cuts shown in Table [I} 
except we require \mui — 50 GeV| < 2.5 GeV and the l^Trei cut is not used. 


sector particles which the Z' can decay into, then the branching ratio to SM states becomes 
~ which can be quite low. The rate for four-lepton hnal states may then be negligibly 

9d 

small. However, there is a much larger rate for only one Z' decay to the SM while the other 
decays to the dark sector, giving a resonance plus MET signal. 

The MET distributions for this signal model are compared to the WW background in 
hgure for two different choices of parameters. The differing cut-offs for the signal distri¬ 
butions motivate multiple signal regions with different MET cuts to maximize sensitivity 
over all of parameter space. As above, using MG simulation we can estimate the reach 
of our proposed search in terms of the total cross-section times branching ratio for the 
(j) —)■ Z'Z' —>■ i'^WxX process, as a function of the </> and Z' masses. We show the results in 
hgure 1^ for two different MET cuts, a 100 GeV cut as considered thus far (Fig. 9a) and 50 


GeV cut (Fig. [9b|). 

In this model the production cross-section of the dark Higgs is determined completely by 
the mixing angle 9h between the dark and SM Higgs states. We can therefore convert the 
above bounds on total cross-section times branching ratio for pp —)■ £~^i~xx to bounds on 
sin^ 9h X BR(Z' —>■ Wi~) (with the approximation BR(0 —)■ i~^i~xx) ~ 2 x BR(Z' —>■ £+£“)). 


The resulting bounds are shown in hgure 10, again for two diherent choices of MET cut. 

















mz' (GeV) 
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Reach in (crx BR)pp^ /+ 100 GeV MET cut 


(Dark Higgs Model) 



W70 (GeV) 



Reach in (crx BR)pp^/+ /-^^(fb), 50 GeV MET cut 
(Dark Higgs Model) 


150 200 250 300 

(GeV) 


(a) 


(b) 


FIG. 9: Left: Projected 95% exclusion reach in the plane of the dark Higgs and Z' masses 
for the dark Higgs model, from an 8 TeV search as in figure]^ with a 100 GeV MET cut. 
In the red region, the decay (j) —)■ Z'Z' is kinematically forbidden. Right: Projected reach 
for the same search with the MET cut reduced to 50 GeV. 


Values of order few times 10 ^ as would be probed by this search could be achieved e.g. 
with sin^ 6^/1 ~ 0.1 (as allowed by Higgs coupling fits jSU [85]) and e ~ O.lgo- 


D. Twin Higgs 

Beyond providing dark matter candidates, specific dark sector models can alleviate the 
hierarchy problem of the SM, if they include states which couple to the Higgs so as to 
cancel the SM Higgs quadratic divergences. This paradigm of “neutral naturalness,” which 
includes models such as Twin Higgs [M] (more generally “Orbifold Higgs” [H3IHH]) and folded 
supersymmetry [H9] , can explain the smallness of the electroweak scale without introducing 
new light fields with SM color charge, which thus far have not been observed. In the Twin 
Higgs model | 86 |, for example, the SM Higgs fields are Goldstone bosons of a global S'f/(4), 
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Reach in Sin^[0;,]xBR2'^/+ 100 GeV MET cut 


(Dark Higgs Model) 



W70 (GeV) 


Reach in Sin^[0;,]xBR2'^/+ 50 GeV MET cut 


(Dark Higgs Model) 



150 200 250 300 
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(a) (b) 

FIG. 10: Left: Projected 95% exclusion reach on the quantity sin^ 6^/1 x BR(Z' —)■ £+£“) for 
the dark Higgs model, from an 8 TeV search as in figure with a 100 GeV MET cut. 
Right: Projected reach for the same search with the MET cut reduced to 50 GeV. 


which is explicitly broken to a gauged SU ( 2 )^ x SU ( 2 )s, where the SU(2)a sector is identified 
with the Standard Model SU{2) sector. If the SU{2)b gauge sector is an exact Z 2 copy of 
the SU{2)a sector, then quantum corrections to the quadratic terms in the scalar potential 
respect the full SU (4) global symmetry and therefore give no contributions to the potential 
of the SM Higgs. Appropriately controlled breaking of the Z 2 symmetry can give realistic 
phenomenology without reintroducing Higgs mass divergences at a problematic level. 

Twin Higgs theories therefore require a dark sector with a multitude of states. The exact 
spectrum of this sector can vary depending on the sources of Z 2 symmetry breaking. The 
dark sector can be very difficult to probe at colliders, as it couples to the SM mainly through 
the Higgs boson and mostly gives rise to invisible final states. However, our generic argument 
that dark sector models offer the possibility of a dilepton resonance + MET signature applies 
in this scenario as well, offering a potential new probe. A kinetic mixing term between the 
U ( 1 ) fields of the two sectors respects the Z 2 (and could easily be generated at high scales). 
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and depending on the spectrnm of the model conld give a signihcant rate for dark sector 
gange bosons to decay to the SM. 

As a specihc example, we will work within the original Twin Higgs model [ 86 ]. To discnss 
the phenomenology however we mnst specify the sonrces of Z 2 breaking. Since the top 
Ynkawa conpling and the SU{2) x U{1) gange conplings give the largest contribntions to the 
Higgs qnadratic divergence, the corresponding conplings in the dark sector mnst respect the 
Z 2 symmetry almost exactly. However, we can break the Z 2 for helds that are weakly conpled 
to the Higgs withont npsetting natnralness (as in e.g the “Fraternal Twin Higgs” [20] )• This 
is necessary to achieve realistic phenomenology, as having many light helds in the dark 
sector snch as the twin electrons, nentrinos etc. is in tension with constraints from Big 
Bang Nncleosynthesis (BBN) as well as cosmic microwave backgronnd (CMB) observations, 
if the two sectors reach thermal eqnilibrinm (see e.g. [91]). One way this can be alleviated 
is if there is only one generation of fermions in the dark sector, corresponding to the third 
SM generation. (A complete generation ensnres anomaly cancellation.) Fnrthermore, the 
Ynkawa conplings of the lighter fermions, the “dark bottom” and “dark tan”, can be larger 
than in the SM withont contribnting too mnch to the Higgs mass. A mass for the dark 
sector photon also helps rednce tension with BBN/CMB bonnds. This can be accomplished 
if we introdnce an explicit mass for the dark sector U{1) (hypercharge) held, denoted B'] 
this soft breaking will not affect cancellation of qnadratic divergences. In the mass basis, 
we will denote the heavier and lighter nentral dark gange bosons as Z' and A' respectively; 
in the limit mB> = 0 these become the exact Z 2 partners of the SM Z boson and photon. 

The most accessible dark sector state at colliders is the dark Higgs, which mixes with 
the SM Higgs as in the model of the previons section. Unlike in the SM, the dark Higgs 
has a “mixed” decay to a Z' and A' when 7 ^ 0, as the mass basis and Higgs conpling 

basis no longer coincide. (For small m^/, the rate for decay to two A'’s is mnch lower.) 
The snbseqnent decays of the Z' and A' then depend on the fermionic spectrnm of the dark 
sector. One possibility is that the lighter A' boson is kinematically constrained to only 
decay to the SM throngh kinetic mixing, while the heavier Z' can decay into dark fermions 
(generically denoted by y). In this scenario dark Higgs prodnction leads to the distinctive 
resonance + MET hnal state, 0 —)■ Z'A', A' —)■ Z' —)■ yy. 

The relevant physical parameters to describe the prodnction and decay rates for this 
process are the dark Higgs mass, the mixing angle 9^ between the two Higgs helds, the 


tha' (GeV) 
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Reach in /- Projected sensitivity to Twin Higgs model 

for Twin Higgs model (m7’=300 GeV) mz'=300 GeV, Sin[0/,]=O.3 



m^{GeV) ff 20 (GeV) 

(a) (b) 


FIG. 11: Left: Projected 95% exclusion reach on sin^ x BR(0 —)■ Z'A' —)■ XX^'^^~) for 
the minimal Twin Higgs model, projected for a search in 100 fb~^ of 14 TeV data using a 
MET cut of 120 GeV. The Z' mass is hxed to 300 GeV. Right: Expected discovery 
significance for the process 0 —)■ Z'A', A' —)■ , Z' —)■ xx fo fhe same search, with mz' 

set to 300 GeV and sinOh set to 0.3. At each point in the {mfj),mA') space, the parameters 
of the model are uniquely determined, so the branching ratios can be computed. 

Z' mass, and the A' mass. These four parameters plus the SM Higgs mass and vev are 
determined by the Lagrangian parameter and the hve parameters of the (renormalizable) 
SM plus dark sector Higgs potential. The dark Higgs branching ratio to Z'A' can also be 
computed from these inputs, so the full rate for 0 —)■ Z'A', A' —>■ Z' —>■ yy is predicted 

given the physical masses and the Higgs mixing. (We approximate the Z' branching ratio 
to dark fermions as unity.) Typical branching ratios for the 0 to this hnal state are of order 
10-3. 

® For a Twin Higgs mode to be completely natural, there must be an explanation for the smallness of the 
scale / where the global 51/(4) (and gauged SU[2)b) is broken, in terms of e.g. strong dynamics or 
supersymmetry (as in [921194] ) near the scale /. Such a UV completion may constrain the parameters of 
the potential or introduce additional non-renormalizable terms. Not being committed to a particular UV 
model, we ignore such effects. 
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The large masses required in this model and the small branching ratio to the desired final 
state limit the discovery potential at the 8 TeV LHC. With sufficient luminosity however, 14 
TeV LHC data can begin to probe this process. As in the previous sections, we simulate the 
signal and the WW, ti and ZV backgrounds at 14 TeV in the presence of pileup, now with 
40 pileup events expected per crossing [95]. We impose the same cuts as in the 8 TeV search, 
except for an increased MET cut of 120 GeV and a modified pileup subtraction algorithm 
Appendix [a|) . Figure 


see 


11 


shows the discovery reach in 100 fb ^ of 14 TeV data, in the 


plane of the dark Higgs and A' mass, for irtz' = 300 GeV. Sizable regions of parameter 
space can be probed at the LHG through this technique, with possible improvements from 
refinements of the search strategy. 

Our analysis here is based on what is essentially still a toy model, as the minimal Twin 
Higgs requires UV completion in order to be truly natural (or for its fine-tuning to even 
be defined). However, our results do suggest that it may be possible to discover the dark 
sector states of Twin Higgs theories by searching for dilepton resonances in association with 
missing energy. Indeed, if Nature stabilizes the electroweak scale through a dark sector, 
then processes of this type are one of very few ways to directly probe the spectrum of the 
theory. 


IV. CONCLUSION 

In general, non-minimality in the dark matter sector opens up new discovery prospects at 
the LHC. In this work, we identified a distinctive collider signature that generically emerges 
when the dark sector contains a U{1) gauge boson. We have argued on general effective field 
theory grounds that such models will often give rise to final states with a dilepton resonance 
from the Z' in association with missing energy. The example simplified models we have 
discussed illustrate the various mechanisms to produce the Z' (in decays or radiation off of 
other dark sector particles), the possibilities for its decay (100% to the Standard Model, or 
mostly to the dark sector), and the different kinematic regimes for the signal (such as highly 
boosted dark sector production through a heavy mediator, or low boost from a compressed 
decay). In much of model space the Z' is completely hidden under SM backgrounds in 
standard resonance searches, but can be discovered through the proposed resonance -|- MET 
analysis. 
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A search for dilepton resonances plus missing energy would be a straightforward extension 
of existing LHC analyses in the dilepton + MET channel. (Indeed, published data [SO] can 
already place bounds, as in Figure [^) We have seen that dark sector models can realize 
a range of resonance masses and boosts. An inclusive analysis which could achieve near- 
optimal discovery reach for this broad range of models would involve two simple extensions 
to existing LHC dilepton -|- MET search selections: 

• Scanning over as wide a range of resonance masses as possible, with the dilepton mass 
window of the signal selection chosen to optimize sensitivity to a narrow resonance 

• Implementing multiple signal regions corresponding to different cuts on MET and 
dilepton px, in particular a signal region with a hard enough cut to almost completely 
eliminate the SM background (which would still retain sensitivity to some signals) 

Such a strategy should saturate the LHC’s potential to discover the Z' plus MET hnal 
state for essentially any dark sector model. Indeed, in forthcoming work other authors [Ml 
l97] also propose a search of this type, motivated by other models. 

Further explorations of non-minimal dark sectors could lead to other interesting signa¬ 
tures. For example, one may consider decays of dark scalars (directly) to the Standard Model 
through mixing with the Higgs. By direct analogy to the discussion we have presented, there 
may exist dark sector processes that could produce these scalars in association with other, 
invisible particles. A search for bb or 77 resonances in association with MET may be able 
to reveal such processes, even if the mixing with the Higgs is too small for the new scalar 
to be observable in direct 2 —)■ 1 production. Note that the Higgs portal and kinetic mixing 
operator along with the neutrino portal HLN exhaust the possibilities for renormalizable 
coupling of a dark sector to the SM. 

If new physics does not couple to the Standard Model gauge bosons, there are a limited 
number of possibilities for probing it at colliders. While identifying resonance provides a 
very powerful tool to pick out signals from SM background, current resonance searches at 
the LHC do not utilize all possible, motivated signatures of possible new physics. Searches 
in hnal states such as resonances plus missing energy can help realize the full discovery 
potential of the LHC. 
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Appendix A: Monte Carlo Simulation 


To derive projected limits on new physics models, we performed Monte Carlo simulation 
of both signal events and Standard Model backgrounds. We generated event samples at 
the partonic level using MadGraph v5 |98]. For the background samples, we used the 
Standard Model UFO model provided with MadGraph. New physics UFO models were 
created using Feynrules v2.0 IHH]. Parton showering and hadronization were performed using 
Pythia v6 [100] . We used Delphes v3 |ini] to simulate detector effects including the addition 
of pile-up (consisting of minimum bias events generated in Pythia), and to reconstruct jets 
(using Fast Jet |102] ). leptons (including isolation cuts), and missing transverse momentum. 
To subtract the neutral pile-up components we used the jet area method |1031 1104] as 
implemented in Delphes (with the “active” area algorithm, and using kx jets with radius 
0.6). To reduce the effect of charged pile-up particles, ATLAS and CMS employ different 
strategies; we implemented Delphes modules to reproduce both of these techniques. In the 
ATLAS selection of ref. |50], a jet is classihed as originating from pile-up if none of the tracks 
associated with the jet come from the primary vertex within the z-vertex resolution. CMS 
in |105] dehnes the parameter (3* for each jet: 




Pt 

tracks Pt 


(Al) 


where '^i^puPp is the pp sum of the tracks associated with the jet that do not originate 
from the primary vertex, while pp is the pp sum of all tracks associated with the 

jet (from any vertex). Pileup jets are characterized by high values of (3*. 

When simulating events at 8 TeV, we ran Delphes with the ATLAS detector card with 
modified charged pileup subtraction as discussed above, and with an average of 21 pileup 


events per crossing. For the 14 TeV projection considered in section HID we simulate an 
















23 



m„ (GeV) 


80 . 





METrd (GeV) 


(a) mu 


(b) ^T^rel 


FIG. 12: Comparison between ATLAS and our simulated background for LHC 8 TeV with 

luminosity of 20 fb“^. 


average of 40 pileup events per crossing. With the increased pileup we hnd that signal vs. 
background discrimination is improved by using a variant of the CMS strategy for charged 
pileup subtraction, identifying jets as pileup if (3* > 0.8 (rather than 0.2 as in the analysis 
of |105j L For 14 TeV simulation therefore we ran Delphes with the CMS card, modified to 
implement this pileup subtraction technique. 

After subtracting pileup tracks and calorimeter energy, we implemented jet finding (for 
jet vetos) and lepton isolation in Delphes to match the object reconstruction in the ATLAS 
search [5^: jets are reconstructed using the anti-fc-r algorithm with Ai? = 0.4, muons are 
required to have the scalar sum of pt of tracks above 400 MeV within AR < 0.3 to be less 
than 16% of the muon pt, and electrons have the same isolation cut plus the requirement 
that the sum of in calorimeter clusters within AR < .3 is less than 18% of the electron 


Pt- 

The main SM backgrounds for our analyses are W^W~ and tt. ZZ and ZW processes 
also contribute (and dominate) for mu near the Z mass. In the signal regions of interest 
to us, the background contribution from processes with fake leptons and/or MET is small, 
as demonstrated in [50]. We computed background event rates by applying selection effi¬ 
ciencies from our MC simulation to the total cross-sections for each process, obtained at 
NLO for diboson processes from [inOj and at NNLL for tt from |in7j (without including any 


dependence of the iF-factors on the experimental cuts). Figure 12 shows a comparison be- 
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tween the results of our MC and the ATLAS simulation for certain background distributions 
from [HU] . 

To obtain bounds on signal models from the published ATLAS results (as in hgure [Sa]) , we 
derive 95% CL bounds from a Bayesian likelihood analysis on the binned data of Figure 3a 
of [50] , with signal distributions as predicted by our MC and a flat prior for the signal cross- 
section times branching ratio. To obtain the projected bounds from the more optimized 
searches we propose, we apply the cuts of Table to simulated backgrounds and signal, with 
modihed J^Trei^ Pt,u' and mu' cuts as discussed in section]^ From the resulting values for 
the expected signal efficiencies and background count rates we can derive the expectation 
value of the 95% C.L. bound on the signal rate one would derive assuming Standard Model 
backgrounds only (as in hgures [^ [^, or the expected signihcance level of the excess in 
this channel if a particular signal model is realized (as in hgure [IT|) . 
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